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Abstract 

A new Lagrangian formalism for self-consistent collective neutrino-plasma 
interactions is presented in which each neutrino species is described as a 
classical ideal fluid. The neutrino-plasma fluid equations are derived from a 
covariant relativistic variational principle in which finite-temperature effects 
are retained. This new formalism is then used to investigate the generation of 
magnetic fields and the production of magnetic helicity as a result of collective 
neutrino-plasma interactions. 
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I. INTRODUCTION 



Photons, neutrinos and plasmas are ubiquitous in the universe |T|J^. During the early 
universe, it is expected that photons and neutrinos interacted quite strongly with hot pri- 
mordial plasmas [Q. Although photons and neutrinos decoupled from plasmas relatively 
early after the big bang there are still conditions today where neutrino-plasma in- 

teractions might be important. For example, during a supernova explosion [^-|^, intense 
neutrino fluxes are generated as result of the gravitational collapse of the stellar core. It is 
generally believed that the outgoing neutrino flux needs to transfer energy and momentum 
to the surrounding plasma in order to produce the observed explosion. 

The self-consistent collective interaction between photons and plasmas is traditionally 
treated classically (i.e., without quantum- mechanical effects), where plasma particles are 
either treated within a fluid or a kinetic picture, while photons are described in terms of an 
electromagnetic field. For a self-consistent treatment of collective electromagnetic-plasma 
interactions (see Ref. for example), one considers both the influence of electromagnetic 
fields on plasma dynamics and the generation of electromagnetic fields by plasma currents. 
The interaction between photons and neutrinos, on the other hand, requires a full quantum- 
mechanical treatment and has been the subject of recent interest 0. 

Neutrino-plasma interactions involve charged and neutral currents associated with the 
weak force PJTO| (through the exchange of and bosons, respectively). The collective 
interactions studied here apply to the intense neutrino fluxes. Discrete (i.e., coUisional) 
neutrino-plasma interactions, on the other hand, involve scattering of individual particles; 
such discrete neutrino-plasma particle effects will be omitted in the present work. 

The purpose of the present work is to investigate the self-consistent collective interaction 
between neutrinos and plasmas in the presence of electromagnetic fields. The inclusion 
of electromagnetic effects is a departure from conventional hydrodynamic models used in 
investigating neutrino interactions with astrophysical plasmas 0. Here, we investigate the 
collective processes 

EM ^ a ^ V (1.1) 

and 

1/ ^ a -> EM. (1.2) 

In the first process, the neutrino [u) dynamics is influenced by an electromagnetic field 
{EM) with a plasma (a) background acting as an intermediary, even though neutrinos are 
chargeless particles. In the second process, electromagnetic fields are generated as a result 
of plasma currents produced by neutrino ponderomotive effects. The problem of magnetic- 
field generation associated with self-consistent collective neutrino-plasma interactions is thus 
investigated here within the context of the process ( |1.2|) . 

A. Notation 

In the present paper, the Latin subscript s refers to different components of the neutrino- 
plasma fluid: the subscript s = v refers to neutrinos while the subscript s = a refers to 
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components of the plasma other than photons and neutrinos. To avoid confusion, we use 
the Greek letters a, /5, ■ ■ ■ for Lorentz indices rather than traditional /U, z/, ■ ■ ■; for example, the 
flux four- vector is J" = Nu"", with proper density (Lorentz scalar) and normalized four- 
velocity u"" = (m°,u). In certain cases, objects with Lorentz indices may not be covariant; 
for instance, the fluid velocity = is not covariant and the number density in a 

given frame n = Nu^ is not a Lorentz scalar. The symbols in bold face are three-vectors 
while those in Sans Serif are four- dimensional tensors (such as F for the electromagnetic field 
strength -F^^). The dot ■ describes the contraction of a Lorentz index or an inner product 
of two three-vectors if in bold face. Here, we employ the metric gap = diag(l, —1, —1, — 1) 
and, hence, a ■ b = aPlP — a • b. 



B. Neutrino Descriptions for Collective Neutrino-Plasma Interactions 

To study collective neutrino-plasma interactions, neutrinos can either be described in 
terms of Dirac spinor fields |]9|-[13|, Klein-Gordon scalar fields classical non-relavistic 

fluids [0, or relativistic quasi-particles |r^|TB[. In all these descriptions, the interaction 
between neutrinos (of type z/) and plasma particles (of species a) is described in terms of an 
effective weak-interaction charge In general, G^y has the following property ||1 1|] : 



Gfjf G (ji/ G (ji/ Gfju^ (■^•'^) 

where a [a) denotes a matter (anti-matter) species and u (z7) denotes a neutrino (anti- 
neutrino) species. The effective charge G^iy depends on the Fermi weak-interaction constant 
(~ 9 X 10~^^ eV cm^), the Weinberg angle dw {sin^ 6w ~ 0.23 IT^), and the species a 
and u. For example, for neutrinos interacting with unpolarized electrons (e), protons (p) 



and neutrons (n), one finds ||TT 



C = V2Gf IS^eSuu. + [h - 2Q^ sin^ 9w) ] , (1-4) 



where /o- is the weak isotopic spin for particle species a (Jg = J„ = —1/2 and Ip = 1/2) 
and Qcr = Q'o-/e is the normalized electric charge. Here, the first term in (|1.4|) is due to 
charged weak currents (and thus applies only to electrons and electron- neutrinos) , while the 
remaining terms are due to neutral weak currents (and thus apply to all species). 

To assist us in investigating self-consistent collective neutrino-plasma interactions in the 
present work, all neutrino and particle species are treated as ideal classical fluids. For this 
purpose, we proceed with the classical fluid limit for plasma-particles in the Dirac description 
expressed in terms of the correspondence 

^,(77c)^. J: ^ in^,Ja), (1.5) 

where ip^j is the Dirac spinor field for particle species a (with 7" denoting Dirac matrices) 
while n„ and J„ = n^v^/c are the particle density and (normalized) particle flux for each 
plasma-fluid species a in the lab reference frame, respectively. In this limit, the propagation 
of a neutrino test-particle of type z/ in a background plasma is determined by the effective 
potential 
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K(x,v,t) ^J2Go 



[1.6) 



where (x, v) denote the neutrino's position and velocity. We note that neutrino propagation 
in matter is a topic at the heart of the problem of neutrino oscillations in matter [pO| -p^ 
and the solar neutrino problem |23|. Although the term J„-y/c is a relativistic correction 
to ria in ( |1.6| ), we keep it for the following reason. For a primordial plasma with a single 



;i-7) 



family of particles (s = a) and anti-particles (s = a), we find from (|1.3|) 

Es=a,W GsuUs = 



and thus the effective neutrino potential ( |1.6| ) becomes Vi, = G^u {Jw — J a) • v/c, for each 
(cr, a) -family. Hence, keeping this relativistic correction is necessary for the description of 
collective neutrino interactions with a primordial plasma P^. The model presented here 
therefore retains all relativistic effects associated with the neutrino and plasma fluids. 

For a self-consistent description of collective neutrino-plasma interactions in which neu- 
trino ponderomotive effects on the background medium are included, we now use a similar 
classical-fluid correspondence for the neutrinos. The propagation of a plasma test-particle 
of species a (with electric charge g^) in a background medium composed of a neutrino fluid 
of type u and an electromagnetic field is determined by the potential 



v;(x,v,t) 



0(x, t) + Y^ G„u n^i^, t) 



A(x, t)+Y^ G^^ J^{yi, t) 



p (1-8) 



where riy and J,^ = n^jVy/c are the neutrino density and (normalized) neutrino flux in the 
lab reference frame, respectively, and A are the electromagnetic potentials, and (x, v) 
denote the plasma-particle's position and velocity. It is interesting to note how the right 
side of (pTSp links the electrostatic scalar potential and the neutrino density n^,, on the one 
hand, and the magnetic vector potential A and the neutrino flux vector J,^, on the other 
hand. We will henceforth refer to the approximation whereby Jo- and Ji, are omitted in 
( |1.6|) and (|1.8| ) as the t(;eaA;-electrostatic (or non-relativistic) approximation. 

Although we assume that each neutrino flavor has a finite mass, this assumption is 
not crucial to the development of our model; see Section || for a discussion of neutrino- 
fluid dynamics for arbitrary neutrino masses. Furthermore we shall ignore all quan- 

in 



I.e., we as- 



tum mechanical effects, including effects due to strong magnetic fields 
sume B / Bqm = hfle/meC^ -C 1, where fie = eB/nieC is the electron gyrofrequency and 
Bqm ~ 4 X 10^'^ G). Hence, although magnetic fields appear explicitly in our model, they 
are not considered strong enough to modify the form of the interaction potentials (|1.6| ) and 



C. Magnetic-Field Generation due to Neutrino-Plasma Interactions 



An important application of the process ( |1.2| ) involves the prospect of generating mag- 
netic fields in an unmagnetized plasma as a result of collective neutrino-plasma interactions. 
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This application may be of importance in investigating magnet ogenesis in the early universe 
(e.g., see Ref. 0]). A similar process of magnetic-field generation occurs in laser-plasma 
interactions whereby an intense laser pulse propagating in a nonuniform plasma generates 
a quasi-static magnetic field. This process was first studied theoretically p6|-p^ and was 
recently confirmed experimentally p9| . 

The generation of magnetic fields by collective neutrino-plasma interactions was first 
contemplated in the non-relativistic (weak-electrostatic) limit by Shukla et al. POBTI]. The 



covariant (relativistic) Lagrangian approach introduced by Brizard and Wurtele |jT5|, how- 
ever, revealed the presence of additional ponderomotive terms missing from previous analysis 
I^,|30|j31| . These additional ponderomotive terms involve the time derivative of the neutrino 



flux dtJu and the curl of the neutrino flux V X Jy (henceforth referred to as the neutrino- 
flux vorticity), which are shown here to lead to significantly different predictions regarding 
neutrino- induced magnetic- field generation. In fact, we show that magnetic- field generation 
due to neutrino-plasma interactions is not possible without these new terms. 



D. Organization 

The remainder of this paper is organized as follows. In Section |T|, the Lagrangian 
formalism for ideal fluids is introduced. In Section |T|, a variational principle for collec- 
tive neutrino-plasma interactions in the presence of an electromagnetic field is presented. 
This Lagrangian formalism is fully relativistic and covariant and can thus be generalized 



to include general relativistic effects (e.g., see Refs. [^,0]). In Section the nonlinear 



neutrino-plasma fluid equations and the Maxwell equations for the electromagnetic field are 



derived. Through the Noether method |34-36|, an exact energy-momentum conservation 



law is also derived and the process of energy-momentum transfer from the neutrinos to 
the electromagnetic field and the plasma is discussed. In Section 0, magnetic-field gen- 
eration, magnetic-helicity production and magnetic equilibrium involving neutrino-plasma 
interactions are investigated. Here, we find that neutrino-flux vorticity (V X J^) plays a 
fundamental role in all three processes. We summarize our work in Section and discuss 
future work. 



II. LAGRANGIAN DENSITY FOR A FREE IDEAL FLUID 

The present Section is dedicated to the derivation of a suitable Lagrangian density for 
a free ideal fluid from an existing single-particle Lagrangian for a free particle of arbitrary 
mass (including zero). The difficulty with dealing with the case of free neutrinos as particles 
is that their mass may be zero. Since the relativistic Lagrangian L for a free single particle 



of mass m is 



L = -mc^- ^ -mc(^^) , (2.1) 

it is not obvious how to handle the limiting case of zero mass. This difficulty is resolved in 
as follows (see also Ref. 0). 
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A. Single- Particle Lagrangian 



Consider the primitive Lagrangian 

Lp = p-±-Ei = -Pacv'^ (2.2) 

for a particle of arbitrary rest-mass m (including zero), where (x,p) are coordinates in the 
eight-dimensional phase space in which the particle moves and = (c, v) = cf Although 
the particle's space-time location = (ct,x) is arbitrary, its four-momentum p° = (i?/c, p) 
is not since the particle's physical motion is constrained to occur on the mass shell 

= m^cl (2.3) 

Here, = 7^ " is the normalized four-velocity and 7 = (1 -|- |up)^/^ is the relativistic factor. 

Since the mass constraint ( |2.3| ) cannot be derived from the primitive Lagrangian ( ^.2[ ), 
we explicitly introduce it by means of a Lagrange multiplier: 



Lp = - p^cv'^ - ^ [m^c" - paP'' ) , (2.4) 



where A ^ is the Lagrangian multiplier and the factor 1/2 is added for convenience. Since 
the Lagrangian (|2.4|) is independent of pa, the Euler-Lagrange equation for pa yields 



f^ = -co'+^^.0, (2.5) 

dpa A 

from which we obtain 

p" = At;"/c. (2.6) 

Using the mass constraint ( |2.3| ) and the identity v ■ v = 7"^, the relation ( p.6| ) yields 

A = 7mc^ (2.7) 

i.e., A is the energy of a single particle of mass m. 

If we now substitute ( |2.6| ) into the primitive Lagrangian (|2.4|) (i.e., by constraining the 
physical motion to take place on the mass shell), we find the physical Lagrangian 



TTl^ A 

L(^;;A) = Lp{x;p = Xv/c;X) = - — . (2.8) 



This Lagrangian now depends only on v" and A (for a free particle, there is no space-time 
dependence in the Lagrangian). The Euler-Lagrange equation for A now yields 

' ' ^ = 0, (2.9) 
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d\ 2 V A2 7 

which gives ( p.7|) . Substituting of ( p.7| ) into ( |2.8| ) yields the standard Lagrangian ( p.l|) . 
For a massless particle, on the other hand, the condition (0|) yields 
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7"' = v^v'' = 0, (2.10) 

which states that massless particles travel at the speed of light. Here, A is still the massless 
particle's energy since ( ^.61 ) gives = A/c. For a massless particle, the single-particle 
Lagrangian is therefore simply given by the last term in (|2.8|), i.e.. 



L{v-\) = -^v^v'^. (2.11) 



This Lagrangian appears in the bosonic part of the Lagrangian for a spinning particle p8 
The Lagrange multiplier A^^ corresponds to the "einbein" which describes the square-root 
metric e = y/g along the world line in a particular gauge where the world line is parameter- 
ized by time. 



B. Lagrangian density for a Free Ideal Fluid 

We now discuss the passage from the finite-dimensional single-particle Lagrangian for- 
malism based on ( p.lj) to an infinite-dimensional fluid Lagrangian formalism. To obtain a 
Lagrangian density for a fluid composed of such particles, we multiply ( p.lj) by the reference- 
frame density n, noting that the proper density is = n7~^. The Lagrangian for a cold 
ideal fluid is therefore 



= -mc^N = -mc^n^/v^ = -mc^^J J'^Ja, (2.12) 

where J" = nf" = (ip'y'^il'/c) is the flux four- vector with a suitable ensemble average (■■■). 
The Lorentz invariance is manifest in the last expression. 

Another contribution to the Lagrangian density of an ideal fluid is the term — Ne{N, S) 
associated with the internal energy density of the fluid in its rest frame, where the internal 
energy e(A^, S) is a function of the proper fluid density and its entropy S (a Lorentz 
scalar). By combining these two terms, the Lagrangian density for a free relativistic fluid is 
therefore written as 

Co = -N[mc'^ + e{N,S)] = -Ne{N,S), (2.13) 

where the total internal energy 

e{N,S) = mc^ + e(A^, ^) (2.14) 

includes the particle's rest energy. 

As discussed above, the single-particle Lagrangian for a free massless particle is given 



as (|2.11| ). The Lagrangian density for a cold ideal fluid composed of massless neutrinos is 



therefore given as 



Co = - Y Ji. ■ = - ^^^^ ■ f i., (2.15) 



where A'^, is a Lorentz-scalar Lagrange multiplier field. The last expression is equivalent 
upon changing the variable Aj^ = n^A^. 
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III. CONSTRAINED VARIATIONAL PRINCIPLE 



The self-consistent nonlinear neutrino-plasma fluid equations presented in this paper are 
derived from the variational principle: 
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J d'x C (A^ F'^^; Ns, <, Ss) = 0, (3.1) 



where in addition to its dependence on the electromagnetic four-potential A" and the Fara- 
day tensor F"^, the Lagrangian density C depends on the proper density A^^ = ns%^, the 
normalized fluid four- velocity = (7<j,Us), and the proper internal energy (per particle) 
Eg for each fluid species s (here, s = a denotes a plasma-fluid species and s = u denotes a 
neutrino- fluid species) . 

The proper internal energy es{Ns, Ss) includes the particle's rest energy [see Eq. ( p.l3D ] 



and depends on the proper density Ng and the entropy 5*^ (a Lorentz scalar). The first law 
of thermodynamics |l39|- ^] is written as 

dEs = TsdS, - PsdN;\ (3.2) 

where Tg is the proper temperature and ps is the scalar pressure for fluid species s. In what 
follows we use the chemical potential for each fluid species s: 

fi, = d{EsN,)/dNs = Es + Ps/N,, (3.3) 

which represents the total energy required to create a particle of species s and inject it in 
a fluid sample composed of particles of the same species. Associated with the definition for 



the chemical potential (3^), we also use the identity 

= T,9"5, + N^'d'^p,. (3.4) 

Note that the independent fluid variables for each fluid species are Ns, u" and Sg although 
other combinations are possible . 

The Lagrangian formulation for the nonlinear interaction between neutrino and plasma 
fluids in the presence of an electromagnetic field is expressed in terms of the Lagrangian 
density 

C = - NsSs - Y.J<r ■ \q.A + Y.GauJu\ + T^F:F, (3.5) 

s=a,u a \ u y IGtT 

where F : F = F'^^ F^^- The first term in (|3.5| ) denotes the total internal energy density of 
fluid s. The second term denotes the standard coupling between a charged (plasma) fluid 
and an electromagnetic field. The third term denotes the coupling between the neutrino- 
fluid species v and the plasma-fluid species o . Note that the second and third terms can be 
written as X^o- ^o-K-, where the single-particle velocity v in ( |1.8|) is replaced with the fluid 
velocity Vg.. The fourth term is the familiar electromagnetic field Lagrangian. 
In the variational principle (|3.1|) , the variation is explicitly written as 

iC . SA.- - ff:- + + Su..— + ^S.— ). (3.6) 
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where 5Fa/3 = daSA^ — d/sSAa so that the second term in (|3l6|) can also be written as 
+2dSA : d£/df. In constrast to other variational principles |]32| , ^ , the Eulerian variations 
SNg, Sug and 6Ss in ( |3.6| ) are not arbitrary but are instead constrained. 

To obtain the correct Eulerian variation, recall that the variation of the fluid motion is 
an infinitesimal displacement of the fluid elements. With a fluid element s described by the 
four-coordinate x", the normalized velocity four- vector is given by 



dr 



dr dr 



-1/2 



dr 



dx'^ 
dr 



(3.7) 



where r parametrizes the world line of the fluid element. Under the infinitesimal displace- 
ment Xg—>-Xg+ [with = {S^g,6$,g)], the apparent variation at a position following 
a fluid element along its worldline is 



dr 



dr 



dr 



dr 



d6^s dxg 
dr dr 



with Us ■ d = \dxs/dT\~'^ djdr and 



(3.8) 



(3.9) 



is a symmetric projection tensor |^0| (i.e., 
space-time location is therefore given by pc 



Us 



0). The Eulerian variation at a fixed 



5u:{x) = hf{us-d)5isp-{5is-d)u', 



(3.10) 



It is easy to check that this variation preserves m"Mq, = 1. 

The variation of the proper density A^^ can be obtained by the requirement that the 
quantity 



Ns 



dr dr 



-1/2 



d'^x 



(3.11) 



should be kept intact (i.e., mass is conserved). The factor in the bracket is the induced 
metric along the world line. This requirement fixes the variation at a position following a 
fluid element along its worldline as —Ns[{d ■ 5^s) — Us/3{us ■ d)6C,'^] = —Ns [hs'^daS^sis], and 
hence the Eulerian variation is given by 



SNs 



i5^s ■ d)Ns -Nshs-. dS^s. 



(3.12) 



It is straightforward to check that the above variations Eqs. ( |3.10| , p.l2| ) are consistent with 
the conservation law 



daJ^ = 



(3.13) 



of the flux four- vector J° = NgU'^. It is useful to know its variation which can be easily 
calculated using Eqs. ( |3.1U| , |3.12| ): 
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5j: = dp{j^5c--m^), 



(3.14) 



where the conservation law ( ^.13 ) has been used. 

Finally, the non-dissipative flow conserves entropy along the world line, 

{us ■ d)Ss = 0. 

To be consistent with the variation Eq. (|3.10| ), we find 

6S, = -{6^s ■ d)Ss. 



(3.15) 



(3.16) 



The expressions ( 3.10 , 3.12 , 3.16| ) give the correct relativistic generalizations of the (non- 
relativistic) constrained Eulerian variations ; see Appendix ^ for a geometric interpre- 
tation of Eqs. (|3.10| , |3.12| , p.l4| , |3.16|) . An alternative variational principle would introduce 
d ■ Js = = Us ■ dSs explicitly in the Lagrangian density by means of Lagrange multipliers 
|3l. 



IV. SELF-CONSISTENT NONLINEAR NEUTRINO-PLASMA FLUID 

EQUATIONS 

We now proceed with the variational derivation of the dynamical equations for self- 
consistent neutrino-plasma fiuid interactions. In deriving these equations, we use the ther- 
modynamic relations ( |3l^ )-( |3.4| ) as well as the continuity and entropy equations ( |3.13| , p.l5| ) 
for each fiuid species s. 



By re-arranging terms in the variational equation (|3.6| ) so as to isolate the variation 
four-vectors 6^s and SA, we find 



sc = d-j - 



dsC + d - I Us ■ - Ns-^ hs 
\ dus oNs 



where dsC = dNg {dC/dNg) + dus ■ (dC/dus) + dSs (dC/dSg), and the Noether four-density 
JT" is expressed in terms of 6^s and 6 A as 

■H.-iV.§^.)-^6.2f .M. (4.2) 

When performing the variational principle ( p.l|) , with SC given by ([4.1| ), we only consider 
variations S^s and 6 A which vanish on the integration boundary. Hence, the Noether density 
J' in ( |4.1| ) does not contribute to the dynamical equations. 

A. Plasma-Fluid Momentum Equation 

First, we derive the relativistic plasma- fiuid four- momentum equation. Upon variation 
with respect to 6^a in (p.lD, we obtain 
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dL ^ dC dC 



+ d 



Ur. 



dC 



- N^- — h. 



(4.3) 



Substitition of appropriate derivatives of the Lagrangian density C and using the constraint 
equations (|3.13| , |3.15|) and the thermodynamic relations ( p.2|) -( p^ , this equation becomes 
the relativistic plasma-fluid four-momentum (covariant) equation 



where 



(4.4) 



(4.5) 



is an anti-symmetric tensor which represents the influence of the neutrino background 
medium g. This tensor satisfies the Maxwell-like equation d^M^^ + d'^M^P + a^M^" = 
and its divergence is d^M^^ = □ , where O = d ■ d and the continuity equation d ■ Jy = Q 
for the neutrino fluid was used. 

Separating the space and time components in (^^) (i.e., using the 3 + 1 notation), the 
spatial components of the plasma-fluid four-momentum equation (|4.4| ) yield 



(9t + v,-V) (/i^7<xvjc2) = -n-^Vp, + (^E + ^ xb) + 



(4.6) 



where fo- is the neutrino-induced ponderomotive force (averaged over neutrino species) on 
the plasma-fluid species a, defined as 



f(T = ^ Go 



Vn^ + ^ ^ j + ^ X V X 



(4.7) 



The neutrino-induced ponderomotive force fo- is composed of three terms: an electrostatic- 
like term Vriy, an inductive-like term dtJu-, and a magnetic-like term V X Jj,. This termi- 
nology is obviously motivated by the similarities with the electromagnetic force on a charged 
particle. In previous work by Silva et al. [|T^, only the electrostatic-like term is retained 
in the neutrino- induced ponderomotive force, i.e., the neutrino particle flux Jy is discarded 
under the assumption of isotropic neutrino and plasma fluids. 



B. Neutrino-Fluid Momentum Equation 



Next, we derive the relativistic neutrino-fluid four- momentum equation; the limiting case 
of zero neutrino masses is treated below ( |4.12| ). Upon variation with respect to in (^.ip, 
we obtain 







dNy 



dC 



+ duy 



dC 



dC 



dUy ' " dSy 



dC 

duy 



•h. 



Nu- — h, 

dNy 



(4.8) 



Substitution of derivatives of C and using the thermodynamic relations ( p.2|) -( p^ , this 
equation becomes the relativistic neutrino-fluid four-momentum equation 
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u,. 



(4.9) 



where 



m5 = - d^J^ 



(4.10) 



is another anti-symmetric tensor which represents the influence of the background medium. 
This tensor satisfies the Maxwell-hke equation d^M^^ + d^Mj^P + df^M^" = and its diver- 
gence is daM^f^ = □ J^, where the continuity equation 9- J^- = for the plasma fluid was used. 
In (|4.9| ), we note that the neutrino fluid is thus under the influence of an electromagnetic-hke 
force induced by nonuniform plasma flows. We also note that the symmetry between the 
ponderomotive forces ( [4.5|) and ( ^4.10|) is a result of the symmetry of the neutrino-plasma 
interaction term {JZaY^v Ga-uJa ■ Jv) in the Lagrangian density (^.51). 

Using the 3 + 1 notation, the spatial components of neutrino-fluid four-momentum equa- 



tion (|]9D yield 



(9t + V^- V) i^Huly^y/c 



(4.11) 



where f,^ is the plasma-induced ponderomotive force (averaged over plasma-particle species) 
on the neutrino-fluid species z/, defined as 



EG. 



(4.12) 



The plasma-induced ponderomotive force on the neutrino fluid is composed of three 
terms: an electrostatic- like term Vn^, an inductive- like term dtJa-, and a magnetic-like 
term V x J„. 

We now discuss the case of a cold ideal fluid composed of massless neutrinos. Variation 
of the neutrino part of the Lagrangian density 

^ ^ G(jv J(j ' Jy 

with respect to b^y yields 

bCy = -UyXy 5Vy ■ Vy " ^ G ^y J „ ' 5 J y , 



(4.13) 



where we used the constraint Vy ■ Vy = 0. Using 6Jy = d ■ {Jy5^y — 5^yJy) and 5vy ■ Jy 
d ■ (JyS^y) ■ Vy, the variation equation (|4.13| ) becomes 



6Cy = d - J + UyS^y " 



Jy ■ d{XyVy) " ^ G ^y M„ " Vy 



(4.14) 



where the tensor Mo- is defined in (^4.101) and the Noether density is 



J = 5^y 



g Gau Jy ■ Ja — Jy {^yVy + ^ Gay J a^ 



(4.15) 
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^From ( [4.14| ) the variational principle / 5Cy d^x = yields the cold neutrino fluid equation 

v,-d{Kv,) = Gau M^-v,. (4.16) 

In the cold-fluid limit, on the other hand, ( [4.9| ) yields u^, ■ d{%^\yUy) = J^a M„ ■ u^, 
where is the neutrino energy. By substituting u,y = ji^Viy into this expression, we readily 
check that ( [4.9[ ) and ( |4.1(j| ) are identical in the massless-neutrino cold-fluid limit and that 



]9D can in fact be used to describe neutrino-fluid dynamics with arbitrary neutrino mass. 

C. Maxwell Equations 

The remaining equations are obtained from the variational principle ( p.lj) upon variations 
with respect to the four-potential One thus obtains 

Substitution of derivatives of C, this equation becomes the Maxwell equation 

d-F = AnJ^q^J^. (4.18) 



Using the 3 + 1 notation, we recover one half of the familiar Maxwell equations from ( [4.181 ). 
The other half is expressed in terms of the Faraday tensor alone as 

QPp'^P + d'^F^P + d^FP"" = 0, (4.19) 
which, using the 3 + 1 notation, yields V • B = and V X E + c'^dtQ = 0. 

D. Energy-Momentum Conservation Laws 

Since the dynamical equations ( |4.3| ), ( [4.8| ) and ( |4.17] ) are true for arbitrary variations 
{S^a-, S^u) and 6 A (subject to boundary conditions), the variational equation (^?T| ) becomes 

6C = d-J, (4.20) 

which we henceforth refer to as the Noether equation. We now discuss Noether symmetries 
of the Lagrangian density ( |3.5D based on the Noether equation (|4.20| ). 



For this purpose, we consider inflnitesimal translations ^ x" + Sx" generated by the 
inflnitesimal displacement four-vector fleld 6x. Under this transformation, the Lagrangian 
density C changes by 

6C = -d-{6xC). (4.21) 

Next, we introduce the following explicit expressions for (5^o-, 6C,y) and 6 A in terms of the 
infinitesimal generating four-vector 6x: 
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6A 



h, ■ Sx 



F-6x - d{A-6x) 



(4.22) 



where the symmetric tensor is defined in ( |3.9| ). (These expressions are given geometric 
interpretations in Appendix]^.) 

Substituting ([4.22|) in the Noether density ([4 .21), we find 



J 



^ dC , 

Us t: 



5x + 2diA-5x) ■ — , (4.23) 



where we have used the identity ■ = in writing the second and third terms. Making 
use of the Maxwell equation ( |4.18| ), the last term in ( [4.23| ) can be re-arranged as 



dC 

2d{A-6x)-^ = d- 



HA-Sx)^ 



(4.24) 



We now note that the expression ioi d ■ J in ( ^.20| ) is invariant under the transformation 
J ^ J + d - K^ where K is an antisymmetric tensor (for which (9^^ K'^^ = 0) which vanishes 
on the integration boundary in ( ^.l] ). Since the first term on the right side of ( [4. 241 ) is 
such a term, it can be transformed away and the final expression for the Noether density is 
therefore 



J 



2 ^ . F - — A + V 



' dC , 
Us ^z— ■ h, 

OUs 



Ns- — h, 

ON, 



■ 6x. 



(4.25) 



Substituting ( [4.21| ) into ( |4.2CI| ), the Noether equation becomes d ■ {J' + Sx C) =0. We 
define the symmetric energy-momentum tensor T from the expression 



J + Sx C = - J ■ 5x, 
where, using ( 4.25 ), the energy- momentum tensor T is explicitly given as 

dC 



T 



dA 



A 



E 



' dC , 

Us 7^ h, 

OUs 



For a constant translation Sx, the Noether equation ( [4.20|) then becomes 

= a-T, 

where using the Lagrangian density ( p.5|) and its derivatives in ( [4.27| ), we find 



(4.26) 



(4.27) 



(4.28) 



1 / „„, g""^ 



EE 



47r I " 4 ■ 



Vsg 



a(3 



(4.29) 



This energy-momentum tensor contains the usual terms associated with an electromagnetic 
field and a free relativistic ideal fiuid [^9| -^. It also contains the energy- momentum terms 
associated with collective neutrino-plasma interactions (third set of terms). 
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The energy-momentum transfer between the electromagnetic-plasma background and 
the neutrinos can now be investigated. Such a process is relevant to supernova explosions, 
for example, where approximately 1% of the neutrino energy needs to be transferred to the 
surrounding plasma. First, we define the electromagnetic-plasma (EMP) energy-momentum 
tensor: 

Temp = ^(f-F-|F:F^+^ {N^fia u„u„ - Pag) = Tem + Tp, (4.30) 



and, using the exact energy-momentum conservation law ([4.28|) as well as the dynamical 
equations (gj), (|]|) and ( ^A8|) , we find 



d ■ Temp = E (^E M,j ■ J^. (4.31) 

This equation describes how energy and momentum are transferred from the neutrinos to 
the electromagnetic field and the background plasma. Note how the transfer of energy- 
momentum between an electromagnetic-plasma and neutrinos is very much like the transfer 
of energy between a plasma (P) and an electromagnetic field (i.e., d ■ Tp = J^a '?o-F ■ J^.) in 
the absence of neutrinos. 



We note that in addition to energy and momentum, wave action [T^ can be transferred 
between the neutrinos and the electromagnetic-plasma background. In this case, electro- 
magnetic waves and/or plasma waves can be excited by resonant three- wave processes. 

V. MAGNETIC FIELD GENERATION AND HELICITY PRODUCTION BY 
COLLECTIVE NEUTRINO-PLASMA INTERACTIONS 

An important application of the process of energy-momentum transfer associated with 
collective electromagnetic-plasma-neutrino interactions is the possibility of generating mag- 
netic fields in an unmagnetized plasma as a result of collective neutrino-plasma interactions. 
Such process might be relevant to the problem of magnetogenesis and the production of mag- 
netic helicity in the early universe [^^-^. A similar process of magnetic-field generation 
has been observed in laser-plasma interactions p6|-p9[]. 

According to our neutrino-plasma fluid model [based on ([4.4|) , (^^ , and (^4.18|) ], the 



strength of the magnetic field generated by neutrino-plasma interactions scales as the first 
power in the Fermi weak-interaction constant Gf- In what follows, we thus refer to magnetic 
fields generated by classical plasma processes (e.g., the Biermann-battery effect and the 
nonlinear dynamo effect) as zeroth-order fields while those generated by collective neutrino- 
plasma interactions as first-order fields. Second-order fields, for example, might be produced 
by processes such a.s a' ^ u ^ a ^ EM, where the first plasma-particle species [a') need 
not be charged (e.g., neutrons). 

In this Section, we investigate the role played by collective neutrino-plasma interactions 
in generating magnetic fields and magnetic helicity as well as magnetic equilibrium. 



15 



A. Magnetic- field Generation 



An equation describing magnetic-field generation resulting from collective neutrino- 
plasma interactions is derived as follows. We begin with Faraday's law 



OB 

'dt 



-cV X E, 



(5.1) 



where for a given plasma-particle species a [using ( |4.7| )], the electric field E is expressed as 



E = -(F. 



-xB, 



where fo- is the neutrino-induced ponderomotive force given by (^l7| ) and 

+ V, • VP, + n-^ Vp,, 



dt 



(5.2) 



(5.3) 



with Po- = (/Uo-/c^)7o-Vo- the generalized momentum for plasma-fluid species a. 

Since the electric field E is common to all charged-particle species, we multiply ( p.2| ) on 
both sides by and sum over all charged-particle species present in the plasma. Defining 
So- 1^ = Q^i the electric field E is then given as 



E 



E 



Substituting explicit expressions for F, and f,, we obtain 



X B. 



(5.4) 



E 



V — 



dtii, - V, X V X + vx<, + v(7;'r,) 



E 



7" cQ^' 



xB, (5.5) 



where 7^ ^T, is the temperature in the lab reference frame and 

Flo- = Pff ~l" X/I/ G 0-1, J yj C 

X,cr — ^ av ~l~ To f^a To ^(jSa 

Eq. ( ^.5|) can then be substituted for the electric field into Faraday's law ( ^.1|) to give 
^ = E || [ V (to^'^o) X - V X (^iP, - V, X V X P.) 



(5.6) 



+ E^Vx(v.xB) - EE 

o V V a 



V X {dtJy - VaXV XJu) 



(5.7) 



The first collection of terms (linear in q^) on the right side of ( |5.7| ) includes the so-called 
Biermann-battery term (Vn^^ X VT^-) ||27|j28|j48[| while the second term (proportional to 
ql) represents the nonlinear dynamo effect. These classical (zeroth-order) terms have been 
known to play important roles in the generation of magnetic fields during laser-plasma 
interactions [p^ -|29[] as well as the evolution of cosmic and galactic magnetic fields 
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The last collection of terms (proportional to QaGo-u) in ( |5.7|) are associated with collec- 
tive neutrino-plasma interactions and are completely new. Here, the neutrino-flux vorticity 
(V X J^) plays a fundamental role in generating first-order magnetic fields; such terms are 
completely missing from previous works ||30| , |3l| . 



According to { f)-7\j , the electrostatic part of the neutrino-induced ponderomotive force 
( [4.7|) does not play any role in generating magnetic fields. Indeed, for each neutrino-fluid 
species u, we have V X [{J2s QsGsu) Vrij^] = 0, independent of the plasma-fiuid composition. 
The neutrino-induced ponderomotive force on plasma particles of species a actually given 
in p!l| , |5n| , pT| is —n~^ {J2s' Gs'^ns') V?7.y = f^^^; this expression improperly involves a sum of 
plasma-particle species {J2s') instead of the sum over neutrino species it appears in 

( ^^ . Shukla et al. then go on to develop a model for magnetic-field generation based 
on the fact that V X f^'^^ 7^ for a plasma with multiple particles species. Since the sum 
over plasma-particle species (I]s') appearing in f^^^ is inappropriate, however, the conclusion 
drawn by Shukla et al. that magnetic fields can be generated in a plasma composed of 
neutrons (cr = n) and electrons (cr = e) by terms such as V(n„/ne) X Vriy is incorrect |^ . 

For a primordial plasma, we note that the Biermann-battery term could be small unless 
the terms V(7~^To-) X VSa and V(7^^Tg^) X V are in opposite directions whereas the 
nonlinear dynamo requires net plasma flow. Using the identities ( |1 . 7] ) , on the other hand, 
we note that particles (a) and anti-particles (a) of the same family (a, a) contribute equally 
to the generation of first-order magnetic fields in a primordial plasma since 



T,s=CT,a Qs Gsu^s = Qa G {^a + V^) ^ 



> . 



(5.8) 



This remark is especially relevant to the problem of magnetogenesis in the early universe. 
Conversely, we note from ( |5.7| ) that a time- dependent magnetic field automatically generates 
neutrino- flux vorticity V X J^- Hence, the usual assumption that the neutrino distribution 



is isotropic fl^] appears to be inconsistent with first-order magnetic-field generation by first- 
order collective neutrino-plasma interactions. 



B. Magnetic Helicity Production 

Another quantity intimately associated with magnetic-field generation is the generation 
of magnetic helicity 

H = I A-Brf^x, (5.9) 

where V is the three-dimensional volume which encloses the magnetic field lines; to ensure 
that this definition of magnetic helicity be gauge invariant, we require that B-fi = 0, where n 
is a unit vector normal to the surface dV . Magnetic helicity is a measure of knottedness (or 
flux linkage) in the magnetic field ^^; hence a uniform magnetic field (or more generally 



a magnetic field which has a global representation in terms of Euler potentials a and (5 
as B = Va X V/3) has zero helicity. The production of magnetic helicity is therefore an 
indication that the spatial structure (and topology) of the magnetic field is becoming more 
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complex. It is expected that this feature in turn plays a fundamental role in the formation 
of large-scale structure in the universe . 

The time evolution of the magnetic helicity (|5.9|) leads to the equation 



dH r f 

— = -2c E-Bd^x - c (0B + ExA) - nd^x, (5.10) 
dt Jv JdV 

where integration by parts was performed in obtaining the surface term. Taking the integra- 
tion volume V arbitrarily large (or requiring that E be parallel to n in addition to B • n = 0), 
we find that the surface term vanishes and we are left only with the first term in ( |5.1CI| ). If 
we now substitute ( ^3|) into ( p.lO|) , we obtain 



dt V Jv 



dtU„ - V,, X V X + Vx. + v(7;'T„ 



d^x, (5.11) 



where Ho- and Xo- are defined in (|5.(j|). Since the term B • Vxo- can be written as an exact 
divergence, it does not contribute to the production of magnetic helicity. Furthermore, 
since temperature gradients along the magnetic field, B-V(7~^To-), vanish in the absence 
of dissipative effects the last term in (|5.11| ) drops out. Hence, magnetic helicity production 
is governed by the equation 

f --i:|?/,B.(W,-v,xVxnj. (5.12) 

This equation states that helicity production can occur in the presence of (zeroth-order) 
non-trivial flows [^| and/or (first-order) nonuniform neutrino flux. 



It has been pointed out that magnetic helicity plays an important role in allowing energy 
to be transferred from small to large scales by a process called inverse cascade. Thus 
neutrino- flux vorticity leads to the generation of small-scale magnetic fields, first, and then to 
the production of magnetic helicity. The production of magnetic helicity, on the other hand, 
converts the small-scale magnetic fields to large-scale magnetic fields which are expected to 
play a fundamental role in the problem of structure formation in the early universe. The 
magnetic helicity production described by (|5.12| ) involves a multi-species fluid picture. A 
more standard description is based on the magnetohydrodynamic (MHD) equations in which 
plasma flows are averaged over particle species. Future work will proceed by deriving ideal 
neutrino-MHD equations. 



C. Magnetic Equilibrium in a Magnetized Plasma and Neutrino Fluid 

When gravitational effects can be ignored, plasmas can be confined by magnetic fields. 
Such an equilibrium is established by balancing the (outward) kinetic pressure gradient with 
the (inward) magnetic pressure gradient. We now investigate how magnetic equilibria are 
modified by the presence of neutrino fluxes. 

The equation for magnetic equilibrium involving magnetic fields associated with neutrino- 
plasma interactions can be obtained by multiplying ( |5.2| ) with qe^n^ and summing over the 
charged-particle species only. In a time-independent equilibrium [d/dt = 0) involving a 
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quasi- neutral plasma (where J2a ^o-^o- = 0), a static magnetic field B and time-independent 
neutrino fiuids, we find the following equilibrium condition 



- X B = V 

c 



J2 (^aV^Pa + I Pa) 



E 



E 



(e g*.. ^) X V X J. 



(5.13) 



where J = (c/in) V X B = J^aQaJa is the current density fiowing in a time-independent 
magnetized plasma. The first term on the right side of ( p.l3| ) represents the classical term 
associated with equilibrium in a magnetized plasma. The second and third terms denote 
first-order neutrino-plasma contributions to magnetic-field equilibrium. 

Shukla et al. ||3^ derived a similar equilibrium condition with only the electrostatic- 
like term present on the right side ( [5.131 ). For a primordial plasma, using ( |1.7| ), we note 
that the neutrino-induced electrostatic-like term once again vanishes from the magnetic- 
field generation picture. Hence, whereas the second term in (|5.13 ) vanishes for a primordial 
plasma, the third term on the right side of ( |5.13| ), however, does not. Magnetic equilibrium 
in a primordial neutrino-plasma is thus described by the balance equation 



(5.14) 



where summation over species on the left side of (|5.14| ) involves only particle species, while 
the summation on the right side involves particle and anti-particle species. Once again, 
neutrino-fiux vorticity V X J^, plays a fundamental role in collective neutrino-plasma inter- 
actions in the presence of an electromagnetic field. 



VI. SUMMARY AND FUTURE WORK 

We now summarize our work and discuss future work. The model for collective neutrino- 
plasma interactions presented in this work is based on the nonlinear dissipationless fiuid 
equations (|4.4|), ( |4.9| ) and ( [4.18| ). These equations are derived from a variational principle 
based on the relativistic covariant Lagrangian density (|3.5|) . An exact energy-momentum 
conservation law ( [4.2(^ ) is obtained by Noether method with the energy-momentum tensor for 
self-consistent collective neutrino-plasma interactions in the presence of an electromagnetic 
field is given by ( [4.29| ). New ponderomotive forces acting on the plasma- neutrino fiuids, 
which are absent from previous works ||l^j3yj31|, are given by ( [4.5| ) and ( |4.1CI| ) [or 



and ([4.12|) , respectively]. In Eqs. (|5.7|) and (|5.13|) , we have demonstrated the crucial role 



played by neutrino-fiux vorticity (V X J^,) in the processes of magnetic-field generation and 
magnetic-helicity production in neutrino-plasma fluids. 

In future work, we plan to further investigate the importance of the new neutrino- 
induced ponderomotive terms associated with neutrino fluxes. For this purpose, it might 
also be useful to derive from ideal neutrino- magnetohydrodynamic equations from ( [4. 41 ), 
( [4.9|) and ( ^4.18|) . Using the new mechanisms for magnetic-fleld generation and magnetic- 
helicity production proposed in ( [5.7D and ( |5.12| ), respectively, we plan to investigate the 
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problem of magnetogenesis in the early universe. As another application, we plan to inves- 
tigate neutrino-plasma three-wave interactions leading to the excitation of various plasma 
waves in unmagnetized and magnetized plasmas; such transfer processes could be important 
during supernova explosions. 
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APPENDIX A: DIFFERENTIAL GEOMETRIC FORMULATION OF 
CONSTRAINED VARIATIONS 



In this Appendix, the geometric interpretation of the constrained variations ( 3.101 , 3.12 



|3.14| , |3.16| ) is given in terms of Lie derivatives along the virtual displacement four-vector 



6C,. Since the variation of a fluid field is only its infinitesimal displacement, all covariant 
quantities are varied by their Lie derivatives with respect to the virtual displacement four- 
vector S^. Here, we use the following definition of the Lie derivative on the /c-form a along 
the four- vector 6C,, denoted L^^a [ p]] : 

hs^a = isc dct + d (i^^ ■ a) . (Al) 

Here, da is a (fc + l)-form while i^^ ■ a is a (fc — l)-form representing the contraction of the 
four- vector 6^ with the fc-form a. By definition, if a = is a scalar field (i.e., a zero-form). 



The constrained variation 6S = — 6^ ■ dS for the entropy S [( p.l6| )] is consistent with its 
geometric interpretation scalar field: 

6S = -Ls^S = -6^-dS, (A2) 

where is^ ■ S = and i^^ ■ dS = {6^ ■ d)S. 

The geometric interpretation of the particle flux J" = Nu" is given as the components of 
the three-form J = {l/?)\)eapK\J°'dx^dx'^dx^ . The constrained variation of the particle-flux 
four-vector is defined as 

5 J = -Ls^J. (A3) 

Since dJ = {d- J) ft with the volume four-form ft = dx^ A dx"^ A dx^ A dx^, and hence dJ = 
due to the continuity equation, we obtain SJ = —d{is^ ■ J), or 

Sr = dpiJf^SC - J'^S^'^), (A4) 



which is Eq. (|3.14| ) itself. /^From this variation, one can easily compute the variations of 



N = y^T^^ and = J'^/N leading QXT^ ) and ( [XTDD , respectively. 
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In Sec. [IV D| , we consider infinitesimal translations — + generated by the 
infinitesimal displacement four-vector 5x. Under this transformation, the Lagrangian den- 
sity C changes \yy 5C = —d- {Sx C). This expression is consistent with the geometric 
interpretation of >C as a density in four- dimensional space, i.e., 

6Cn = -Ls.{Cn), (A5) 

where is the Lie derivative with respect to 6x. Here, using isx ■ d{C f2) = and 

d [isx -{1^^)] = d{C6x-uj) = d- {6x C) n, (A6) 

we easily recover ( [4.21| ). 

Next, the expressions for 5 A is given in ( [4.22| ). Here, the electromagnetic four-potential 
A appears as the the components of the one-form A ■ dx. Thus 

5A-dx = -liSx{A-dx). (A7) 

Since ■ d{A ■ dx) = — (F ■ 5x) ■ dx and d[iix ■ {A ■ dx)] = d{A ■ dx), we easily recover ( |4.22| ) 
for the four-potential A. We note that the expression = h-Sx given in ( [4.22|) is consistent 
with the expressions 6S = — L^^S = — L^S* and 6J ■ lo = — h^^^J ■ uj) = — Lsx{J ■ i-o). 
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